# Multivariate Regression Analysis (Tables 4A, 4B, 4C)

model1 <- lm(support[trust, understanding]~group,data=data)
model2 <- lm(support[trust, understanding]~group+sex+age+education+ethnicity+income,data=data)
model3 <- lm(support[trust, understanding]~group+sex+age+education+ethnicity+income+party,data=data)
model4 <- lm(support[trust, understanding]~group+sex+age+education+ethnicity+income+party+ideology,data=data)
model5 <- lm(support[trust, understanding]~group+sex+age+education+ethnicity+income+party+ideology+military,data=data)
stargazer(model1,model2,model3,model4,model5,type = "html", style = "qje",
          digits = 2,
          ci = TRUE,
          title = "Support[Trust, Understanding] for AI in Different Domains and for Different Purposes", 
          covariate.labels = c("T1 (Cars, Enhance)","T2 (Cars, Substitute)","T3 (Online, Enhance)","TS (Online, Substitute)","T5 (Drones, Enhance)","T6 (Drones, Substitute)","Sex","Age","Education","Race",
                               "Income","Political Party","Political Ideology","Military Service")
          dep.var.labels   = "Support",
          df = FALSE,
          omit.stat = c("rsq", "ser"))